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Abstract 

Brane inflation in superstring theory ends when branes collide, initiating the hot big bang. 
Cosmic superstrings are produced during the brane collision. The cosmic superstrings produced in 
a D3-brane-antibrane inflationary scenario have a spectrum: (p, q) bound states of p fundamental 
(F) strings and q D-strings, where p and q are coprime. By extending the velocity-dependent one- 
scale network evolution equations for abelian Higgs cosmic strings to allow a spectrum of string 
tensions, we construct a coupled (infinite) set of equations for strings that interact through binding 
and self-interactions. We apply this model to a network of (p, q) superstrings. Our numerical 
solutions show that {p, q) networks rapidly approach a stable scaling solution. We also extract the 
relative densities of each string type from our solutions. Typically, only a small number of the 
lowest tension states are populated substantially once scaling is reached. The model we study also 
has an interesting new feature: the energy released in (p, q) string binding is by itself adequate 
to allow the network to reach scaling. This result suggests that the scaling solution is robust. 
To demonstrate that this result is not trivial, we show that choosing a different form for string 
interactions can lead to network frustration. 
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I. INTRODUCTION 


The cosmic microwave background (CMB) data from WMAP strongly supports the 


inflationary universe [3(] as the origin of the hot big bang. Recently, the brane world scenario 
suggested by superstring theory was proposed, where the standard model of the strong and 
electroweak interactions are open string (brane) modes while the graviton and the radions 
(which govern the size/shape of the bulk) are closed string (bulk) modes. Consider a generic 




brane world scenario that closely describes our universe today (i.e., a KKLT-like vacuum 
If we take an extra brane-anti-brane pair in the early universe their brane tensions provide 
the cosmological constant that drives inflation. There is an attractive force between a brane 


ejnflation is taking place. 
7|; in brane inflation. 


and an anti-brane so they tend to move towards each other whi 
Thus, brane inflation is a natural feature of the brane world 5 
the inflaton is an open string mode identified with the interbrane separation. Inflation ends 
when the D3-D3-brane pair collides and annihilates, releasing energy that starts the hot big 
bang. Note that the inflaton field no longer exists after the annihilation of the D3-D3-brane 
pair. Towards the end of inflation, the D3-D3-brane collision produces Dl-branes - or cosmic 
superstrings - but neither monopoles nor domain walls j^. 

We can estimate the cosmic string tension /i using the density perturbation magnitude in 


the CMBR data from COBE 


Q 


In the simplest realistic scenario, namely, the KKLMMT 


D3-D3-brane inflationary scenario [ 3 ], one finds that 


5 X 10"^ > G/i > 4 X 10 


-10 


( 1 ) 


where G is the Newton’s constant. The upper bound comes from WMAP data and other 
data ^1, while the lower values require some fine-tuning in the model. These predic¬ 
tions are not sensitive to different warping schemes because the normalization to COBE is 
always performed after the warp effect is taken into account. In this scenario, the density 
perturbation responsible for structure formation is dominated by the inflaton, with cosmic 
strings playing a secondary role. 

The evolution of networks of cosmic strings is a well studied problem ^ Q]. After 
the initial production of cosmic strings, the strings interact among themselves. When two 
cosmic strings intersect, they reconnect or intercommute. When a cosmic string intersects 
itself, a closed string loop is broken off. Such a loop will oscillate quasi-periodically and 
gradually lose energy via gravitational radiation. Its eventual decay transfers the cosmic 
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string energy to gravitational waves. A higher (lower) string density leads to a higher 
(lower) interaction rate so, not surprisingly, cosmic string networks evolve towards scaling 
solutions. A consequence of scaling is that the physics of simple, abelian Higgs networks is 
essentially dictated by a single parameter, the dimensionless string tension Gfi. This scaling 
feature can be seen by considering the evolution of the number density in a one-scale model, 
where the scaling solution emerges as an attractive hxed point. 

There are many different hybrid inflationary models in which one can construct a variety 
of cosmic string-producing scenarios. What is different in brane inflation is that the string 
network that is produced has a large spectrum of possible string tensions [l4 , [3, • For 

the D3-D3-brane inflationary scenario, one expects a spectrum of {p, q) string bound states 
[3, [l^, where the tension of a particular bound state (p, q) is given by 

Gfi{p,q) = Gp^Jp^g1 + q^. ( 2 ) 


where Ps is the superstring coupling. For the string to be stable, p and q must be coprime, 
i.e., p and q have no common factors greater than 1. Written this way, (1, 0) corresponds to 
the fundamental F-string while (0,1) corresponds to the Dl-brane, or D-string, so p is the 
tension of the (0,1) superstring. We can see the emergence of such a spectrum in several 
ways. Consider the following simple picture: the gauge group at the end of inflation just 
before brane-antibrane annihilation is 


U{l)n X 17(1)5 = U{1)+ X 17(1)_, 

where the open string (complex) tachyon held stretching between the branes couples only 
to 17(1)_. The usual operation of the Higgs mechanism generates abelian vortices following 
the spontaneous symmetry breaking of 17(1)_. These are the D-strings. Since no free 17(1) 
gauge symmetry remains after the annihilation, it is believed that the 17(1)_|_ symmetry 


becomes conhning, yielc 
strings, or F-strings 


ing conhning huxes that may be identihed as fundamental closed 
17j |. The production of D-strings may be estimated via the Kibble 
mechanism. Most of the decay products are expected to be very massive non-relativistic 
closed strings, which are expected to decay to gravitons, standard model particles and other 
light modes. We expect some of the massive closed strings to be extended. These are the 
F-strings. In a cosmological setting, their production is likely, again, to be dictated by the 
Kibble mechanism. The production of D- and F-strings are not independent, so we expect 
some initial spectrum of (p, q) strings to be produced. 
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Since the interactions between (p, g) strings is not simple [18|, one expects that (p, g) 
network evolution might be quite involved. It is not obvious, a priori, that the network can 
even approach scaling. For example, it could oscillate (i.e., the density of any specihc (p, g)- 
type could oscillate indehnitely), approach scaling only asymptotically, or simply frustrate. 
One way to address this problem would be to do a full numerical simulation of a (p, g) 
string network. However, this is a highly non-trivial problem. String network evolution is a 
complex physical process; accurately modeling the build-up of small-scale string structure is 
computationally demanding, even in the context of abelian Higgs models, which have only 
one type of vortex. A radically simpler alternative would be simply to generalize the one- 


scale string network model due to Kibble [l^ to the case of (p, g) string evolution. Recall 
that the scaling of the cosmic string network appears as a stable hxed point in this one-scale 
model. However, previous researchers have found it useful to include more of the network 
physics than the original one-scale model allowed. In particular, the velocity-dependent 
one-scale (VOS) model developed by Martins and Shellard for the abelian Higgs case Q 
provides a very convenient and reliable method for calculating the large-scale quantitative 
properties of string networks in many contexts, including a cosmological setting. This model 
performs exceptionally well when tested against high resolution numerical simulations of 


string networks It allows one to see analytically how scaling emerges, and to calculate 
reliably a small number of macroscopic quantities useful for cosmological applications. We 
take this model as the starting point for our own model building. We recognize that there are 
a number of other analytic approaches to the string evolution j^. Some may characterize 
the details of small-scale stringy structures more accurately, but they also require more 
phenomenological input parameters which can only be obtained from simulations. Since 
there is, as yet, no simulation of a cosmic superstring network, and since we are chiefly 
interested in the overall properties of such a string network, we choose the simplest possible 
“analytic” model that highlights the most important physical effects. 

In this paper we adapt the velocity-dependent one-scale model to describe a multi-tension 
cosmic string network that includes both string self-intersection and string-string binding 
interactions. For a multi-tension network, the string density evolution equation generalizes to 
a set of (inhnitely many) coupled equations. We then specify the string interaction terms in 
our particular multi-tension, (p, g) string model and solve this set of equations numerically. 
Fortunately, we hnd that the (p, g) string network, with stringy interactions turned on. 
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rapidly approaches scaling. This fast convergence and rapid decrease of string densities 
with increasing tension allows us to truncate the set of equations at low, computationally 
tractable values of {p,q)- To show that this scaling result is not somehow built into the 
model we have constructed, we demonstrate that the same set of evolution equations with 
a different interaction term can lead to a frustrated network. 

Because of the various approximations we use, our results are limited to overall macro¬ 
scopic network features, which are nonetheless those features that are needed for cosmological 
applications. For instance, we hud for the {p, q) superstring network : 


• The {p, q) string network approaches a scaling network rapidly. The hnal scaling 
solution is independent of the initial densities of the various types of strings. The 
fractional density in strings, for F 7 ^ 0 , is given by 




20/(0.55P + F) ^, = 1.0 
15/(0.53P + F) ^, = 0.5 


( 3 ) 


where P measures the probability of self-interaction and F measures the overall prob¬ 
ability of interaction among different types of strings. For the {p, q) cosmic superstring 
network, we do not know the value of F, though we expect P < F < 1. It is inter¬ 
esting to note that scaling is achieved even if we turn off the string self-interaction, 
i.e., when P = 0. For the abelian Higgs model we have Ft/(i) ~ 20 and P ~ 0.28. 
Using this value for P and taking F = P, we hud that, for = 1, F ~ 46; for F = 1, 
P = 0.28, we have F ~ 17. Thus, the total density of the (p, q) cosmic superstring 
network is comparable to standard cosmic strings. Differentiating between the two 
kinds of networks based on their string densities will require more detailed modeling. 


• The relative number density of each type of string is roughly given by 


~ F(m) 6 < n < 10. (4) 

The fall-off is power-like, not exponential. The rapid convergence of the coupled set 
of equations is brought about by this rapid fall-off. The power law is most accurate 
for high values of p and g; the spectrum tends to be somewhat flatter for the hrst few 
string types. Indeed, we hud that when scaling is reached, the relative numbers of 
( 0 , 1 ), ( 1 , 0 ), and ( 1 , ± 1 ) strings are comparable and far larger than the population of 
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the remaining (p, q) states with p, \q\ > 1. In the case of F = P = 0.28, ps = 1.0, we 
find A^(p,g) oc 

• The adapted multi-tension velocity-dependent one-scale model (MTVOS) that we de¬ 
scribe in l^can be used for many different kinds of multi-tension networks by a simple 
change of the inter-string interactions term. 


Clearly, this analysis can be improved in a variety of ways. We shall comment on some 
of them. However, we are conhdent that the rapid approach to scaling and the fast power 
drop-off in the densities are generic features of cosmic superstring networks. 

In the past, the evolution of cosmic strings in models more complicated than the abelian 
Higgs model have been considered. For example. Pen and Spergel studied the evolution of S'tv 


strings by simulating a network of S 3 and Sg strings 


22j |. The S'tv symmetry enforces identical 


tension and number densities among the N string types. In terms of the model we construct, 
the set of N coupled evolution equations collapses to a single equation. McGraw also modeled 


non-abelian S 3 strings with two different tensions 


The evolution of networks of 


strings connected to monopoles, w’ 


11c 


consider, has been studied in Ref. 


1 have some qualitative similarities to the networks we 
24l |. We believe that a network of (p, q) strings is the hrst 
network type that truly requires a set of coupled equations. The formalism may be adapted 
for other non-trivial string network. 

In Section 2, we adapt the velocity-dependent one-scale model to a model for the evolution 
of comic strings that have a spectrum of tensions. In Section 3, we specialize the model to the 
(p, q) superstring network by dehning our string interaction term. In Section 4 we present 
our numerical results and in Section 5 we briefly discuss some observational implications of 
these networks. 


II. THE MULTI-TENSION VELOCITY-DEPENDENT ONE SCALE MODEL 

Consider a set of different types of cosmic strings {a} with tensions {pa}- Let the number 
of cosmic strings of type a per unit area be Uq. Suppose that all of the cosmic strings may 
be characterized by a single length scale L and a single average velocity v, and that cosmic 
strings of type a can evolve either by interaction mediated loop formation or by binding 
to cosmic strings of other types f3 ^ a. The following model is motivated by the model of 
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Martins and Shellard [^, bnt has been altered somewhat to accommodate the new string 
physics we introdnce. 

We assnme that the length scale evolves via the eqnation 

L = HL + ciV , (5) 


where the loop parameter ci < 1 is a dimensionless factor and H is the Hnbble parameter. 
We take the eqnation of motion for the velocity to have the Martins-Shellard form 


V = {1 -v^) (^-2Hv + ; 


( 6 ) 


where the “momentnm parameter” C 2 is a second constant. This term is the acceleration 
dne to the cnrvatnre of the strings. In the absence of expansion, these two eqnations imply 
7 = (1—= {L/, and with ci = C 2 = 0, bnt with expansion retained, they imply 
-yva^ = constant, L/a = constant, and so, = constant. Thns, the “self-acceleration” 

dne to string cnrvatnre and expansion have opposite effects: self-acceleration increases string 
velocity, whereas expansion dilntes it. This snggests that the two effects can cancel one 
another. We can see how this comes abont by rewriting Eq. © as 


d{HL) 

dt 


H {HL + civ) + —HL = H{HL) 


CiV 

fl + 3w\ 

HL 

V 2 )\ 


( 7 ) 


and combining it with the velocity eqnation; we then End that there is a qnasi-steady solntion 


v = HL 


n + 3w\ 

V 2ci ) 


C2 

2HL ’ 


( 8 ) 


which implies 


HL = 


C2C1 

l-t-3w’ 


V = 


1 Ic2{1 + 3w) 

Cl 


( 9 ) 


In this solntion, both HL and v are constants that differ in the radiation {w = 1/3) and 
matter {w = 0) eras. Clearly, there is no qnasi-steady solntion for w < —1/3; thns, qnasi- 
steady solntions only exist in the radiation and matter eras. We reqnire that in both eras, 
n < 1, a condition that is more restrictive in the radiation era, where it demands that 
C 2 / 2 C 1 <1. In practice, we choose ci and C 2 snch that the scaling valnes oi HL and v 
match the valnes given in j^, where similar constants c and k are chosen to line np with 
fnll network simulations. The translation between our constants and theirs is simple: The 
translation between our constants and theirs is simple: in the radiation era, our C 2 = k, 
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while our ci = (A; + c)/2); in the matter era, C 2 = (3/4)fc and ci = (3/8)(fc + c). For example, 
in the radiation era they hnd HL = 0.1375 and v = 0.655, which for us hxes ci = 0.21 and 

C2 = 0.18. 

Next, add to these equations an equation of energy conservation, at hrst in the absence 
of interactions between strings of different types. Let the cosmic string energy density be 

n/i 

^ \/l — 

where fi is the mass per unit length of a string, and n is the mean string number density. 
In the absence of interactions, assume that 



p = -2Hp{l + 


( 11 ) 


Differentiating the expression, Eq. oa for the string energy density and using Eq. mj and 
Eq. (jni), the equation for h, implies 

n = - ( 2 H + n . (12) 


where we see that the hrst term comes from cosmological expansion. The second term can 
be interpreted in two different ways. For straight strings, it would rehect a net expansion in 
the velocity held orthogonal to the strings. For kinky strings, we should interpret n as the 
characteristic number of intersections per unit area on average for a two dimensional surface 
intersecting the network. As the strings straighten out (the source term in v) the number of 
intersections will fall, at a rate that is ~ nn/L characteristically. We may interpret this to 
mean that as string kinks straighten out, the number of intersections of strings with any two 
dimensional cut through three dimensional space will decrease; it is also the straightening 
of kinks that raises v. Note that when the string velocity approaches its asymptotic value 
V = C2/2HL, the energy equation becomes 


p=- {2H + 


C2V 

~Tj 


V = - 


H 


2 + 


€ 2(1 + 3w) 

2ci 


P , 


(13) 


i.e. it assumes exactly the same form as the equation for the number density of strings. This 
is sensible because the energy per unit length per string /ly —>■ constant asymptotically. 

Next, let us consider what happens when we allow interactions between strings. Recall 
that we assume a single characteristic length scale, L, for all types of strings; however, 
unlike the single-/! case, this length need not be related directly to string density. We 







further assume that the different types of cosmic strings interact by binding, first at a point 
and then zipping up to form a new cosmic string with the same length as the original two, 
which enforces equal lengths for all different kinds of strings. The zipping up takes a time 
L, but let us suppose that L is small enough compared to the Hubble length that we can 
regard the zipping up as instantaneous; HL ~ 0.1 is good enough for our purposes. Let the 
equation for Ua be 


2 ^ ^ 

/3,7 d,7 

(14) 

where the hrst term on the RHS arises from the breaking off of loops from individual undu¬ 
lating strings of type a, the second term arises from breaking off of loops after the collision 
of two strings of type a, and the third term arises from the zipping up of two strings of 
different types that collide and bind. 

As an aside, we note that the term proportional to P could equally well be written as a 
term under the summation, proportional to Poaa- We have pulled this term out of the sum 
to make the contrast between self interaction and interaction between strings of different 
types clear. But a few comments on this term are in order: 


ha + ^Hua =-p- Pn\vL+FvL 

Jj 


• The term proportional to P or Poaa is the usual term that drives networks without 
multiple string types to scaling. 

• By writing P instead of Poaa, we are assuming that the self-interaction rate does not 
depend upon a. 

• If we take Pa /37 ~ ^ li then a self-interaction term of this form drives all string 

species which have a sufficient initial density to the same hnal scaling number density. 

• After we have taken P out of the sum, we either restrict the sum to /? 7 ^ 7 or assume 

Poaa 0 . 


We assume that our constants C 2 , P, and P are identical for cosmic strings of all types. 
We define P as a measure of the overall probability that two strings of different types interact 
at all. We have assumed that the interaction of strings of two different types can only result 
in zipping up, if anything - there are no reconnections and no breaking off of loops directly 
associated with such interactions. 
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In the third term, Pq /37 is the probability of forming a string of type a when strings of 
types (3 and 7 collide, whenever the strings interact at all. The factor 

- (1 + 5^^) 

is introduced so that we do not double count the production of strings of type a when strings 
of different types f3 and 7 collide; i.e. since we do not restrict the sum, symmetry on /?•«-»• 7 
implies we have two identical source terms from /? + 7 —>■ a; the factor 


1 + 


in the loss rate arises because in each a — a collision we lose two long strings. 

By dehning Na = ofua, we can rewrite Eq. m as 

- ^ Pai3jNf3N^{l + Sf3^) — ^ Pi3^aN-yNa{l + 6^a) , 

/ 3,7 / 3,7 

(15) 

and if we dehne conformal time by dr] = vdt/a and introduce the comoving string length, 
L = ai, then we hnd 


Na = 


C 2 N^v PNlvL ^ FvL 


nL 


C2Na 

I 


- PNli + Pi 


- Pa0^NisN-y{l + dp^) — P0^aN-yNa{l + d^a) 
/3,7 /3>7 


(16) 


with a prime denoting differentiation with respect to conformal time. In terms of iV^, we 
hnd that pa = PaNa/ofyl — ^. 

The remaining two equations are those for ^ and v. Substituting L = ia into Eq. © 
gives 

i=— ^ i = £(0) + ci7 . (17) 

a 

When we change the independent variable from f —^ 7 , Eq. © becomes 

= 11 ^ UHav + ^i) . 

To complete the set of equations, we need to hnd 0 ( 7 ); we use 

^ (i(lna) d(lna)d ?7 vd{lna) 
dt drj dt a dr] 


(18) 

(19) 


to get 


da Ha^ 

dr] V 


( 20 ) 
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In the radiation dominated era, which is of greatest interest to us practically, Ha^ is ap¬ 
proximately constant. Thus, in the quasi-steady state, the scale factor grows linearly with 

T]. 

Eq. (IT^ may yield a steady state solution of the form Na = fa!^^ = /q/[^( 0) + Cip]^ 
provided that 


- 2ci/„ 


-{c2U + Pfl) + F 


9 X] PalS'rfldf-f (1 + Sl3'y) 
- A7 


^ ^ Pfi'jaf'yfa (1 T ^"(a) 
0,1 


( 21 ) 


has a nontrivial solution. 

It is instructive to consider what we get when F = 0. In that case, Eq. m has two 
solutions, /q = 0, which is not relevant, and 


/a 


2ci — C 2 

P 


( 22 ) 


which is physically realizable only if 2ci > C 2 . Let us assume that this is so. At sufficiently 
late times, we will therefore find that 


Pa(2ci - C 2 ) 

PcKariy^/l — 


(23) 


where we have assumed that Cip ^ f*(0), and that v relaxes to its asymptotic value. In this 
limit, we find that a ~ Ho?r]/v in the radiation era, and therefore r] ^ va/Ha^ = v/P[a] use 
this in Eq. (Pn|) to find 

^ IJ,aH^{2ci - C 2 ) _ SnGpa ^ 87iGpa{2ci - C 2 ) , . 

Pc\v‘^\/l — a/I — 

Because this version of the network equations assumes common L and v for all string types, 
when P = 0 we expect to find Qa/Pa independent of a, assuming nonzero initial populations. 
(Remember that hlo oc /„ = 0 is also a solution.) Notice that Eq. implies 12^ cc 

(2ci — C 2 )/P, and for small self-interaction probability, this is the expected P~^ scaling. 
Moreover, nonzero P is essential for time-independent hlo, to arise in networks with only 
interactions among strings of the same type. Also, since this is the limit in which our model 
reduces to the abelian Higgs case, we can use prior simulation results to fix the value of our 
parameter P. Numerical studies of the radiation dominated era have T = flcs/ (SvrG/i) ~ 20, 
which for us implies P = 0.28, taking ci = 0.21, C 2 = 0.18, and v = 0.655. Because of this 
we take P = 0.28 as the fiducial value for P in our numerical solutions. 
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It is also instructive to consider what happens if P = F = 0. In this case, Eq. (d has 
an exact solution 


Na{0)m 


cijcx 

N = ... 

" [£( 0 ) + 


iVa(O) 


m 

CiT] 


cilc\_ 


(26) 


(Eq. is a special case of the general solution 


1 

1 pe{Q)' 


Na ~ 

_Na{fi) 2Ci-C2_ 

hc))] 


C2/ci 


+ 


pe 


2ci - C 2 

that can be found for P = 0, and the results of this and the previous paragraph follow from 
appropriate limiting cases of this solution.). From this result, we hnd that 

87rGp„iV„(0)[£(0)]^^/‘=i 


fin 


(26) 


3(Pa)^ c2/ci.yC2/ciy^]^Tr^ 

In the radiation dominated era, when Ha oc a~^, we hnd that Hq cc which either 

rises or falls depending on the sign of 2 — C 2 /C 1 . 

Note that the network equation, Eq. m, may also be used to describe entanglement. 
In that case, instead of setting a = (paHa), as we shall do to describe {p,q) networks, we 
simply let a = Pa- Moreover, we set PaP'y = 5o-(/3+-y); taking C 2 = P = 0, the network 
equations are 


N' 

^ ’ OL 



^a-(l3+'i)Ni3N^{l + 5j3^) — Na -^^(l + ^'ya) 

Pn 7 


(27) 


These equations cannot lead to a scaling solution because there is a conservation law, basi¬ 
cally conservation of energy, that restricts the evolution of the system. Multiply Eq. 
by a and sum over a; the result is 



E aNaN^ 


a, 13 


Fi 

Fi 

0 , 


- 5^(/3 + 7)iV^iV, - aN^N, + 2 /5iV| - 2 aiV; 

- /3^7 


«7^7 


/3 


Y,PN,n,-Y, aN^ 

./37^7 (^7^'y (3 OL 


(28) 
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where the next to last line was obtained by relabelling 7 —/3 in one of the two sums over 
7 ^ 7 . Thus, here we have an example where scaling is not achieved. The network evolves 
toward ever larger values of /i, but its overall comoving energy density does not decline. 

Thus, neither the existence nor the nature of a scaling solution for a particular multi¬ 
tension network is obvious. If each type of string evolves independent of all other types, 
scaling will be achieved eventually for all types present originally , with Hq oc /Xq/P. Turning 
on the interactions between string types will populate the different tensions, and, once 
produced, their self-interactions and energy-losing binding interactions will propel them 
toward scaling solutions. The hnal spectrum of string tensions may be broad or narrow, 
depending on the efficiency with which the reaction terms operate. The reaction terms 
themselves may promote scaling even if there are no self-interactions, but the fact that 
entanglement can be described by a reaction network with particular choices of interaction 
probabilities shows that there are certainly circumstances in which scaling cannot arise solely 
from the reactions among strings of different types. Note, finally, that we have assumed very 
little about the nature of the multi-tension network that is described by these equations 
beyond the assumption that string interactions lead either to loop formation or the formation 
of other kinds of strings through some sort of binding. Thus these equations may easily be 
adapted for any particular multi-tension string network model simply by determining the 
form of Pa /37 for that model; the particular (p, q) network that we consider below is only one 
example of the sort of network these equations can describe. 

III. F- AND D-STRING NETWORK 

To specialize the preceding network to the (p, q) strings of [ 3 ,I 3|5 we define Pa /37 ~ taking 
a = {p,q), P = {k,l), and 7 = {m,n) - and motivate the overall interaction probability, P. 
For this investigation, we make a first and very crude approximation: we assume that the 
probability of two strings of different types interacting is a single, universal constant, rather 
than a function of of a,/ 3,7 or the relative velocity of the strings. By discarding all these 
complexities, we retain only a kinematically determined branching ratio (see both Fig. ^ 
and the discussion in the text below); in future work, we may attempt to retain more of 
the physics contained in P to obtain more realistic results. Before we can write down this 
branching ratio, we shall state the relevant properties of (p, q) strings, since these determine 
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the form of Pap'y- 


• Strings with positive and negative values of p and q are generically allowed; the sign of 
p or q indicates the direction of the string’s charge. Because of a reflection symmetry, 
we can always choose the orientation of the string such that p > 0. For p = 0, g > 0; 
for p > 0, g G Z. 

• Strings with g = 0 are only stable for p = 1 . It is probable that (0, g) strings are 
marginally bound; operationally, we assume that the non-zero momentum transfers in 
the string collisions that accompany binding unbind these states. An interaction that 
formally would create an {N, 0) or (0, N) string thus, in fact, creates N (1, 0) or (0,1) 
strings. 

• We assume that two strings of different types interact with probability F. If two 
strings of different types do interact, there are two possible products, or bound states, 
that they can form; a (p, g) string interacting with a (p', q') string can form either a 
(p + p^ g + gO string or a (p — p', q — q') string, where we always take p > p'. As stated 
above, if either of these product bound states has a resulting p” and q” that are not 
coprime, then what is actually formed is a set of strings with stable, lower p, g values. 

• In agreement with our comments above, we assume Poaa = 0 or, equivalently, restrict 
our summations to /? 7 ^ 7 . 

• For bound states of strings to be stable, p and g must be coprime. If a bound state of a 
string is formed with p and g not coprime, then the new state is, in reality, a collection 
of lower-tension p and g strings that are coprime: i.e., a string which nominally has 
p = Nk, q = Nl is actually a set of N {k, 1) strings. We may view this as the “decay” 
of a {Nk, Nl) string: 

— comes about because the N {k, 1) strings that compose this “state” are BPS with 
respect to each other; that is to say, they have no mutual binding interactions. 
Any possible marginal binding may be ignored since the strings are moving with 
relativistic speeds. 

— has no energy cost - the resulting collection of strings always has lower energy 
than the strings that bound to form them; there is no actual {Nk,Nl) bound 
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state - a collection of stable {k, 1 ) strings is formed immediately following the 
interaction. 

— can untie itself through reconnection events if the collection of {k, 1 ) strings that 
are created in the interaction are tangled or tied up immediately after their 
creation. 



FIG. 1: A schematic view of a string intersection. The intersection angle, 6 , determines whether 
the additive - {p + p',q + q') - or subtractive - {p — p',q — q') ~ binding occurs. 

Whether the interaction of two strings forms a (p+p', q+q') bound state or a (p—p^ q — q') 
bound state is determined by a simple consideration of force balance: if the angle between the 
interacting strings is small enough, then the heavier, additive bound state is formed because 
the two interacting strings’ tensions can balance the tension of the heavier bound state; if 
the angle is greater than some critical angle of force balance, then the lighter, subtractive 
bound state is formed. The critical angle which determines which binding occurs is given 
by 0 

cos emn = ([m - Cn]gs, n) 

\^kl \ |^mn| 

where C is the RR scalar. If we assume a stochastic distribution of string orientations, then 
the strings’ interaction angle should have a flat distribution in cos 6 , that is, that each value 
of the cosine between -1 and 1 should be equally likely. If we assume this, and remember 
that the directionality of the F and D charge must be taken into account - i.e., 6* = 7r/4 
is not equivalent to 6^ = Svr/d - then the probability of forming the additive bound state is 
simply the fraction of cosine-space with 6 less than the angle of force balance; the subtractive 


15 




bound state is formed otherwise. We can write this as: 

p± {k + Cl){m + Cn)g^,+ln 

a/37 2 V ^ V[(^ + Cl^ds + + Cnyg'l + n2]V2 

where indicates a = (p, g) = (A: + m, / + n) and P~ indicates a = {p,q) = {k — m,l — n). 
This form captures the kinematic branching ratio, but as stated leaves out an important 
process: the creation of non-coprime (p, q) strings, which are in reality collections of two or 
more coprime strings. To take this into account, we must slightly modify the way in which 
we insert Pap^y into our equations: we take 

P{Nk,Nl)(p,q){p',q') = ^-P(fc,/)(p,g)(p',g') • 

The inclusion of this process is extremely important. This break-up of non-coprime strings 
is a nonreversible process that is fundamentally dissipative - it helps to keep the average 
tension of the network low both by limiting the pathways by which high-tension bound 
states can be reached and by providing a mechanism through which a single interaction can 
destroy a high-tension bound state and replace it with a collection of low-tension strings. 
Thus, in summary, 

• Two strings of different types a and (3 interact with probability F. 

• When these strings interact, either a subtractive or additive bound state is formed, 
with probabilities P~ and P^ given by Eq. (12911 . The two interacting strings are 
annihilated in the production of the new bound string state. 

• When the bound state 7 is stable, one such string is produced. 

• When the nominal bound state 7 is unstable, it immediately forms N lower-tension 
stable strings, which helps keep the tension dependance of string density spectrum 
steep. 

In all our numerical runs we have assumed C = 0. 

Some possible physical effects that our model neglects include: 

• Velocity dependence of the interaction probabilities due to variation in the binding 
energy of the resulting bound states: the increase binding energy that holds very high 
tension states together is small - the energy gained by binding decreases greatly for 
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high-yU states. Thus we expect that the momentum transferred in even moderate- 
velocity interactions that involve these lightly-bound states may lead them to unbind 
spontaneously. 

• This velocity dependence may lead to an effective cut-off in p, irrespective of the 
interaction dynamics. 

• We might not expect strings of widely varying tensions all to have the same velocities 
and characteristic length scale in reality, contrary to what our model assumes. 

• We have decoupled the evolution of the r.m.s. velocity and length scale of our network 
from the network’s P and F dependent interactions, though one would generally expect 
these interactions to be relevant to determining the network parameters. 


IV. NETWORK RESULTS 

The equations given in m require numerical solution. For all numerical results, we work 
in the radiation-dominated era, assume (for convenience) that the RR scalar, (7 = 0, and 

n 

hx our constants Ci = 0.21 and C 2 = 0.18 to match j20[ (these choices are made so that, at 
scaling, we have HL = 0.1375 and v = 0.655). Furthermore we have done each run twice, 
with two different values of the superstring coupling, Qs = 0.5 and 1.0. We were less certain 
about how to initialize the cosmic string network; cosmic string creation immediately after 
brane inflation is not understood completely. Fortunately, scaling has proven to be quite 
robust to a wide variety initial conditions; for an illustration of this, see Fig. |21 On energetic 
grounds, we believe that, in general, networks will be formed with primarily the lowest-lying 
states populated; thus, for our calculations we chose initially to populate only the (1, 0) 
and (0,1) states, and those with equal number densities. Final scaling results are always 
insensitive to these choices; at worst, very different initial conditions can alter the rate at 
which the network approaches the scaling regime. We similarly hnd that any initial choices 
for network velocity, u, and length scale, L, quickly approach their analytically-predicted 
scaling values. To integrate our equations numerically, another choice we had to make was 
how many (p, q) states to allow. After testing networks of many different sizes, we found that 
our results showed a steep, power law dependence of number density on tension in all cases. 
The relative densities of the low-lying tension states, furthermore, were not changed when 
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more high-tension states were included. Thus we were able to obtain accurate results from 
a relatively small network: for the runs we show here, we have taken p G [0, 5], g G [—5, 5], 
though the way in which we solve the equations allows nominally higher-valued, temporary 
(p, q) states to “form” if the values of p and q are non-coprime, but only if the decay products 
of the unstable (p, q) state are a collection of stable {k, 1) strings with k, |/| < 5. Finally, we 
take the scale factor of the universe a = 1 at network initialization. 



FIG. 2: Comparison among three different sets of network initial conditions, all taking F = 1, 
P = 0. The higher three lines represent the evolution of the overall density in cosmic strings 
(summed over all string states), (ics = f^cs/((8/3)7rGp(op)), with scale factor, a. The lower three 
lines represent the evolution of the density in (0,1), or D-, strings, flis’ ')=b!S’'V((8/3)7rGp(o,i)), 
with a. Our standard initial conditions, equal initial populations of (1,0) and (0,1) strings and 
HL = 1, are shown by the dashed lines. The solid line represents the results from a network run 
with a short initial length scale (10“^ of our usual choice) and with over half of the initial string 
(p, q) states in our network equally populated. Finally, the dotted line shows the results for a very 
large initial population of strings - Qcs ~ 1000 ~ equally spread over half of the tension states 
included in our network, with our usual choice for the initial network length scale. For all the runs 
shown here we have set the superstring coupling, = 1.0. 

An interesting new result from this network model is that string networks with no loop 
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creation - those with P = Poaa = 0 - still exhibit cosmologically acceptable scaling; enough 
energy is lost through string binding and binding-mediated annihilation to keep the comoving 
network number densities constant, regardless of initial conditions, after an initial re¬ 
laxation period following network formation. Because of this, these networks are very robust: 
though there are regions of parameter space where the network never truly reaches scaling, in 
all reasonable cases (where we keep F 7 ^ 0 , P < F) we never hnd cosmologically-disastrous 
solutions where cosmic strings come to dominate the energy density of the Universe. 

There are three regimes of interest for these solutions: 

1. For F S> F : The network will be dominated by F, D, and (1, ±1) strings. Higher 
tension states are present, but maximally suppressed (these networks have the steepest 
spectra). 

2. For F —> 0: All string tensions that are present initially eventually reach the same 
scaling density. If there are a great many string states, this can cause a catastrophe, 
since the formation of loops tends to drive all types of strings to the same value of 
f^cs/ 

3. For F ~ F: The interactions terms will populate the higher [p, g)’s, and the F terms 
will flatten the spectrum somewhat because of its tendency to equally populate all 
levels. The larger F is, the more quickly this happens. In Figs. 0111 El &Ela variety 
of combinations are shown. For larger values of F, the “hnal” scaling state is not an 
exact scaling solution, but one that continues to evolve slowly to late times. Some 
features that appear to generic in this regime include 

• Hcs/((8/3)7rG/r(o,i)) = 60/(F -|- 0.55F) for gg = 1.0, 46/(F -|- 0.53F) for gg = 0.5. 
This formula is only valid for F 7 ^ 0. 

• “scaling ~ 1000 ’ scaling is achieved at T^calmg ~ ^^~%ehea.P "^^^re T^eheat 

is the temperature to which the universe reheats at the end of inflation. Since 
^reheat ~ string scale, ~ 10“^Ms TeV, so scaling is reached 

long before the electroweak phase transition. 

• Scaling results are insensitive to initial conditions unless F = 0 

• Steep hnal spectra: oc /i“", 6 < n < 10 

There are several aspects of these results that require further discussion. 
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A. Scaling 

• The overall properties of the networks are fairly insensitive both to initial conditions 

and to particular parameter choices; i.e., hies never grows fast enough ever to come 
to dominate the universe. 

• However, the hnal state of the string network depends upon the relationship between 
F and P. When P F, the network does not quickly reach a true scaling solution. 
Instead, it continues to evolve to late times (see Figs. ElandlSI). Because of the efficient 
energy loss from both binding and loop formation, the network’s overall density grows 
very slowly ~ e.g., {d log fl/d log a) ~ 0.07 for a ~ 10"^; ~ 0.01 for a ~ 10®, for the case 
of P = F = 0.28; note that {dlogQ/dloga) = 0 dehnes entry into the scaling regime. 
This late growth comes about because of the continuing competition between loop 
formation, which wants equally to populate all string states, and binding interactions, 
which tend to destroy high-tension bound states. This late evolution is not dangerous 
cosmologically. 

• In agreement with our understanding of the late evolution of P ~ F networks, such 
networks tend to develop somewhat flatter hnal spectra, as their high-tension bound 
states tend to be more populated than those in P ^ P networks. Their spectra still 
exhibit very steep power law behavior, however (see Figs. EJandE]). 

• The scale factor at which the network enters the scaling regime is somewhat dependent 
upon initial conditions: networks with more states initially populated tend to take 
slightly longer to reach scaling, though the greater frequency of interactions caused 
by such initial conditions means that these networks tend to be less dense throughout 
their evolution than networks that are formed with only low-lying string states (see 
Fig. 121); networks that begin with much smaller initial L, on the other hand, can 
take a good bit longer to reach scaling. Recall also that the binding interactions of 
high-tension states will very often lead to non-coprime combinations, which leads such 
networks quickly to develop the same kinds of steep spectra that are seen when less 
democratic initial conditions are used. 
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B. Low-F Catastrophe 


• An aspect of superstring networks that has been as-yet unappreciated is that their 
ability to populate arbitrarily high tensions through the formation of bound states 
can lead to a cosmological catastrophe if such states cannot be made to decay. In 
traditional network evolution, which is what our equations reduce to when F —>■ 0, 
even very small initial populations of each possible string state will each eventually 
reach the same hnal scaling density. When this happens, flcs = 

can become huge, and thus disastrous, even if the energy density in each individual 
state is small. 

C. Final Spectra 

• The fact that our numerical solutions have found a very strong dependence of string 
number density on tension - with oc /i”"", and 6 < n < 10 - is another important 
aspect of these networks. If the spectrum were flat, or nearly so, a scenario very 
much like the low-F catastrophe outlined above would ensue: since the effect of many 
string states is additive, and since there are many more possible states at higher 
tensions, such a flat-spectrum network would be ruled out immediately by cosmological 
considerations. 

• Computationally, the steepness of the spectra greatly eases our task. Numerical tests 
showed that the addition of many high-tension states with low number densities 
scarcely affected any of the results. We were thus able to limit ourselves to small 
networks, with pmax = |gmax| = 5. 

• Careful study of Figs. El and El shows that the spectra plotted there are not strict 
power laws. The spectra are, in fact, somewhat flatter for very low tensions, where 
the number of possible states is small. Thus, in all cases there are proportionally more 
F, D, and (1,1) strings than anything else; particularly when F 3> P, these states 
will dominate the cosmic string network. The relative populations of these low-lying 
states are tabulated in Tab. HI 

• The effect of varying the superstring coupling, Qg, is to vary the tension of the F 


21 




gs = 

= 1.0 

gs = 

= 0.5 



M(0,1) 

M(0,1) 

_ 1 M(l.±l) _ 
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F 

P 
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A"(i,±i)/A"(o,i) 

1 

0 

1 

0.769 

3.24 

1.21 

1 

0.14 

1 

0.803 

2.50 

1.19 

1 

0.28 

1 

0.836 

2.19 

1.19 

0.56 

0.14 

1 

0.829 

2.24 

1.19 

0.56 

0.28 

1 

0.887 

1.96 

1.21 


TABLE I: The relative populations of the three lowest-lying tension states, which in all cases 
dominate the networks’ energy density. 

strings relative to the D strings, with this variation propagating up the ladder of 
bound states. Here, we have taken only two representative values of Qg'- 0.5 and 1.0. 
Reducing gg affected the network as our iV oc /x results suggest: those states which 
were previously degenerate (equally populated)- (1, 0) or [(2,1) (2, —1)] relative to 

(0,1) or [(1,2) -|- (1, —2)], for instance - had their degeneracy lifted, with the lighter 
state’s number density increasing. The precise amount of increase was dependent upon 
the values of F and P. Again, see Tab. Q] 

• We note in passing that (unphysical) networks with P = 0 and with the loop term 
oc C 2 removed from Eq. m also go to scaling. 

• In all our networks, the lowest tension states dominate the network energy density. 
We expect this to be a feature of any multi-tension networks that interact via binding, 
even if the spectrum of possible bound states is much more complicated than the one 
we have considered. 


V. OBSERVATIONAL CONSEQUENCES 

Several aspects of (p, q) networks are potentially observationally distinct from regular 
cosmic string networks. The most obvious difference is that these networks feature a spec- 
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FIG. 3: The bottom panel shows the evolution of lies = ^^cs/((8/3)7rG/X(o,i)) for various parameter 
values, taking the string coupling gg = 1.0. The top panel shows the rate of change in the comoving 
number density Nrj'^; in the scaling regime, dlogNrf /dlogg = 0. 


trum of string tensions. We suggest a few possible observational signatures that could allow 
one to distinguish a (p, q) network from a standard, abelian Higgs network: 

• Previous studies of cosmic string lensing probability Q have been based on results 
from standard, abelian Higgs network models. In such networks, P = Vies/(SnGfi) ~ 
20. In principle, our model allows values for P both less than and greater than the 
abelian Higgs value. However, we expect the extra-dimensional nature of superstrings 
to reduce their interaction rates, which leads to higher values of P and Hes- If future 
lensing surveys hnd a rate of cosmic string lensing substantially higher than that 
predicted by a abelian Higgs network, that rate could both be a signature of a cosmic 
superstring network as well as an observational constraint on the parameters of such 
a network. 


If the overall densities in cosmic superstring networks are generally higher than those 
in abelian Higgs models, then observational bounds on cosmic string tension (e.g. 
that depend on overall network properties will need to be reinterpreted. We expect that 
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FIG. 4: The final scaling-era spectra for a variety of parameter combinations, taking the string cou¬ 
pling gs = 1.0, with normalized to unity and the other number densities altered accordingly. 


the net effect will be to tighten such bounds, though how much the bounds will change 
is difficult to predict since observational bounds depend on many aspects of string 
networks (e.g., string substructure, or “wiggliness”), while the detailed properties of 
multi-tension networks have not yet been fully fleshed out. 


The F-shaped junction of two strings in the act of binding is a good signature of 
the non-trivial properties of cosmic superstrings. Such a junction, if present, could 
be detected by cosmic string lensing, or by observation of the Kaiser-Stebbins effect, 
where a temperature difference is seen in the cosmic microwave background radiation 


due to a string-induced Doppler shift 


30|. In the latter case, we would expect to see 


a different temperature in each of the 3 patches of sky. The relativistic motion of 
the binding strings could also be an indicator of a binding event: the cosmic string 
lensing angle is enhanced by a factor of 7 for moving strings (depending on string 
orientation), which is moderate for usual network motions (~ 1.3, in the radiation- 
dominated era). The strings motions near a binding site, however, are very relativistic, 
though over a very small spatial region, and thus would exhibit exaggerated lensing 
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FIG. 5: The bottom panel shows the evolution of lies = ^^cs/((8/3)7rG/X(o,i)) for various parameter 
values, taking the string coupling gg = 0.5. The top panel shows the rate of change in the comoving 
number density Nrj'^; in the scaling regime, dlogNrf /dlogg = 0. 

near the binding site. Random variation among string velocities within the network 
should also lead to the existence of some individual fast-moving strings whose lensing 
will also be enhanced. 

• A recent analysis Q of the direct detectability of cosmic string-generated tempera¬ 
ture anisotropies in the data from the upcoming Planck satellite suggests that rela¬ 
tively high tension or fast moving strings within cosmic superstring networks would be 
marginally within the Planck range of detectability (they estimate that strings with 
Gjji ^ 6 X 10“®, V = 1/v^ would be directly detectable by Planck; for a multi-tension 
network with a fiducial tension Gfi ~ 5 x 10“^, only strings with a combined and 
high-(p, g) tension enhancement of ~ 10 would be seen). 

• Direct observation of more than one cosmic string tension from observational tech¬ 
niques, such as gravity wave bursts Q or gravitational lensing, that are sensitive to 
a particular string’s tension would be a definite prediction of this kind of network. In 
the case of lensing, however, the velocity, orientation, and string substructure depen- 


25 








1 



FIG. 6: The final scaling-era spectra for a variety of parameter combinations, taking the string cou¬ 
pling gs = 0.5, with A^(o,i) normalized to unity and the other number densities altered accordingly. 


dences of the lensing angle may overwhelm this effect for the most probable lensing 
strings (since over 90% of the strings in onr network are F, D, and (1, ±1) strings, 
whose tensions are all of the same approximate magnitnde). For random string ori¬ 
entations, the string lensing angle can vary by as mnch as a factor of six becanse of 
velocity dependent effects, thongh we expect typical velocity dependent variation of 
only a factor of two or so (these variations arise becanse the string lensing angle is 
proportional both to the sine of the orientation angle of the cosmic string relative to 
the line of sight as well as 7(1 + n ■ v), where v is the string velocity and h is a nnit 
vector in the direction between the observer and the string 3l|). Thns, there are two 
ways that lensing measnrements conld indicate the existence of a mnlti-tension net¬ 
work. The most dramatic would be a single very large (> 10) variation between two 
observed lensing angles. Just as compelling, however, would be if a large number of 
lensing measurements were made with a typical variation among events that is greater 
than one would expect to arise, statistically, from random string orientations and 
string velocity directions. In any event, accurate follow-up observations of the Kaiser- 
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Stebbins effect, where the string’s relative velocity enters the equations differently (in 
a cross product rather than a dot product) could perhaps allow us to disentangle to 
some degree the string’s velocity and intrinsic tension. Another possible avenue for 
discriminating between variation due to string orientation and velocity and intrinsic 
string tension would be if a series of lensing events were observed along a single long 
string within a small patch of the sky. Since we expect strings to be curved, it could 
be possible to observe the same cosmic string at several different orientations. This 
could allow us to extract the string’s actual tension through a statistical analysis of 
the events’ lensing angles. 

• Another signature would be a mismatch between a particular string tension measured 
directly ~ from lensing, perhaps - and a Gfi measurement coming from a technique 
like CMB fluctuations or pulsar timing analysis that is only sensitive to the averaged 
network as a whole; however, the effects of string substructure (i.e., string wiggliness), 
which alter string-generated CMB spectra, could mask more subtle expressions of this 
effect. If limits on string substructure (see H) improve, then a direct detection via 
lensing of a string with a tension that is several times larger than what CMB limits 
would lead us to expect for a single-tension network would be a strong indication of 
the existence of a multi-tension network. 

By combining different observations and accumulating sufficient data, one should be able 
to measure a set of properties of the cosmic strings and so distinguish between different 
scenarios. This goal will be easier to reach if the true cosmic string tension is closer to 
today’s observational bound. 
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